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If 7-ray bursts have a cosmological origin, the sources are expected 
to trace the large-scale structure of luminous matter in the universe. 
I use a new likelihood method that compares the counts-in-cells dis- 
1/^ . tribution of 7-ray bursts in the BATSE 3B catalog with that expected 

' from the known large-scale structure of the universe, in order to place 

0^ , a constraint on the distance scale to cosmological bursts. I find, at the 

■ 95% confidence level, that the comoving distance to the "edge" of the 

^ ' burst distribution is greater than 630 h~^ Mpc {z > 0.25), and that 

, the nearest burst is farther than 40 Mpc. The median distance 

' to the nearest burst is 170 Mpc, implying that the total energy 

released in 7-rays during a burst event is of order 3 x lO^'^ h~^ erg. 
None of the bursts that have been observed by BATSE are in nearby 
galaxies, nor is a signature from the Coma cluster or the "Great Wall" 
likely to be seen in the data at present. 
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INTRODUCTION 



The origin of 7-ray bursts is still unknown and is currently the subject of a 
^n , "great debate" in the astronomical community. Do the bursts have a galactic 
origin (|^) or are they cosmological (|^) ? And what is their distance scale? 

In this paper, I do not attempt to answer the first question, but rather, I 
show that if one assumes that 7-ray bursts are cosmological in origin, one 
Q ' can begin to answer the second question and place a constraint on the dis- 
' tance scale to the bursts. This is because cosmological bursts are expected to 
■ trace the large-scale structure of luminous matter in the universe . The 
constraint comes from comparing the expected clustering pattern of bursts on 
the sky — which will depend on their distance scale because of projection 
effects — with that actually observed. The observed angular distribution is in 
fact quite isotropic (||) ; hence, only a lower limit to the distance scale can be 
Ci placed because a sufficiently large distance will always lead to a sufficiently 
isotropic distribution on the sky. 

Here I use a powerful new likelihood method ^ , which I had previously 
developed to analyze repeating of 7-ray bursts in the BATSE IB and 2B 
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catalogs (|6|) , to compare the observed counts-in~cells distribution in the new 
BATSE 3B catalog (|^) with that expected for bursts at cosmological distances. 

Here I will assume for simplicity that fig = 1 and A = 0, and that the 
large-scale structure clustering pattern is constant in comoving coordinates. 
The results are in fact insensitive to these assumptions because of the small 
redshifts that are involved. I follow the usual convention and take h to be the 
Hubble constant in units of 100 km s^^ Mpc^^. 



LIKELIHOOD METHOD 

Let A'ccii be a large number of circular cells, each centered on a random 
position on the sky. Each cell is of fixed solid angle size 17 = 27r(l — cos^iad), 
where 6'rad is the angular radius of the cell. I set the number of cells to 
be such that any part of the sky is covered, on average, by one cell; hence, 
Aceii = 47r/r2. Let Cn to be the number of these cells having N 7~ray bursts in 
them, out of the A'^tot = 1122 in the BATSE 3B catalog, where A^ = 0, 1, 2, ... 
I then define the observed counts-in-cells distribution, Pjq = Cn/NcsIu as 
the probability that a randomly chosen cell of size has A^ bursts in it. The 
counts-in-cells distribution contains information about clustering of 7-ray 
bursts on scales comparable to the angular size ^i-ad of the cell. 

I now define Qjv to be the counts-in-cells distribution that is expected if 
7-ray bursts are cosmological in origin and trace the large-scale structure of 
luminous matter in the universe. This expected distribution depends on only 
one unknown parameter, the effective distance D to 7-ray bursts (which I 
define below), because the angular clustering pattern of bursts on the sky will 
depend by projection on this distance. 

The likelihood C measures how likely it is that the observed counts-in- 
cells distribution Pjv is drawn from the expected distribution Qpf. Since Qn 
depends on the unknown effective distance D to 7-ray bursts, the likelihood 
is really a measure of how likely a given value of D is. I find that (||) : 

log C = Nccw ^ Pn fog Qn + constant . (1) 

N 

Now the cumulative Cmax/ Cnun distribution of 7-ray bursts seen by BATSE 
begins to roll over from a —3/2 power-law for bursts fainter than Cmax/Cmin ~ 
10. Since this is many times above threshold, it suggests that BATSE sees 
most of the source distribution and that this distribution is not spatially 
homogeneous . I define D as the comoving distance beyond which the 
source density drops appreciably. It is not the distance to the very dimmest 
burst in the BATSE catalog, but rather the typical distance to most of the 
dim bursts in the sample; thus, D is the effective distance to the "edge" of 
the source distribution in the BATSE catalog. 

I take the power spectrum which characterizes the large-scale clustering of 
7-ray burst sources to be the same as that determined from a redshift survey 
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of radio galaxies . This power spectrum is characteristic of moderately 
rich environments, and is intermediate between that of ordinary galaxies and 
clusters. Because the exact bias factor relating the clustering of 7-ray burst 
sources to that of luminous matter is unknown, such an intermediate ansatz 
is reasonable. In any case, the resultant distance limit depends only weakly 
on the bias factor (roughly as the square root). 

Knowledge of the power spectrum permits a calculation of the expected 
angular clustering pattern, the expected counts-in-cells distribution Qn, and 
finally the likelihood C [from equation (1)], all as a function of the effective 
distance D to 7-ray bursts (||). I have included the smearing due to finite 
positional errors on the clustering on small scales ( p^ . Indeed, each burst in 
the BATSE catalog is assigned a positional uncertainty 0crr corresponding to 
a 68% confidence that the true burst position is within an angle ^crr to the 
position listed in the catalog. 

I have chosen the cell size ^lad in order to maximize the sensitivity of de- 
tection, or signal-to-noise, given the strength of the signal expected. For a 
sample of 1122 bursts (the total number of bursts in the BATSE 3B catalog) 
with positional smearing of 0^^ — 3.8° (the median value in the 3B catalog), 
the signal-to-noise is maximized when cells of 6'i.ad = 5° are used (P) . 

RESULTS 

Figure 1 shows the likelihood of the BATSE 3B catalog data as a function of 
the effective comoving distance D, calculated using cells of size ^lad = 5°. The 
likelihood is normalized to that expected for an isotropic distribution on the 
sky. At large values of D (the maximum value allowed is D = i?H = 6000 h~'^ 
Mpc, the size of the horizon in a closed universe), the likelihood goes to 
unity, because by projection a sufhciently large distance will always lead to 
an isotropic distribution on the sky. Note also that there is no value of D for 
which the likelihood is greater than 1; thus, the maximum likelihood value for 
D is i?H and the 3B data are consistent with isotropy. 

The solid line in Fig. 1 shows the likelihood for a positional smearing of 
6'orr — 3.8°, corresponding to the median value in the 3B catalog. To illustrate 
the dependence of these results on positional errors, I also show (dashed line) 
the results for a larger positional smearing of 6'orr = 6.6° (with cells of size 
^rad = 9° to maximize signal-to-noise). 

Small values of the effective comoving distance to 7-ray bursts are un- 
likely, according to Fig. 1: I find, at the 95% confidence level, that for the 
3B median positional error of ^orr = 3.8°, D must be greater than 630 ft.^^ 
Mpc, corresponding to a redshift z > 0.25. If the positional errors are larger 
than quoted and are better characterized by 6'orr — 6.6°, these results are 

■^Graziani & Lamb (jl^) compare the 3B positions with those from the IPN network, 
and conclude that the systematic errors are larger than the 1.6° value quoted in the 
3B catalog. Their best-fit model gives a median positional error of 6.6°. 




FIG. 1. The likelihood of the BATSE 3B catalog data as a function of the effective 
comoving distance D to 7-ray bursts, shown with a smearing of ^err = 3.8° and 
6.6°. 



only slightly weakened: At the 95% confidence level, D must be greater than 
500 Mpc, corresponding to a redshift z > 0.19. 

These limits arc not sensitive to earlier assumptions on cosmology and clus- 
tering evolution since these only become important at higher redshifts. They 
are also conservative limits in that a constant median value for the positional 
errors was used, rather than the entire distribution of errors. This is because 
the bright bursts, which ostensibly are nearer to us, are more clustered and 
are responsible for the bulk of the expected signal, but in fact have smaller 
errors than the median value. The faint bursts, which are far away, are hardly 
clustered to begin with (even before smearing), but have errors larger than 
the median value. Hence the expected clustering pattern has been smeared 
more by using a constant median value (this permits a simpler calculation) 
than by smearing using the entire distribution of errors. So the counts-in- 
cells statistic has been weakened somewhat and thus the quoted lower limits 
are in fact conservative. 



CONCLUSIONS 



If 7-ray bursts are cosmological and trace the large-scale structure of lumi- 
nous matter in the universe, and their positional errors are as quoted in the 
3B catalog, then the lack of any angular clustering in the data implies that 
the observed distance to the "edge" of the burst distribution must be farther 
than 630 Mpc. Since there are 1122 bursts in the catalog, an effective 
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limit on the nearest burst to us can be placed by convoluting the likelihood as 
a function of D (Fig. 1) with the nearest neighbor distribution of 1122 bursts 
inside a sphere of radius D. I find that the nearest burst must be farther 
than 40 Mpc at the 95% confidence level, and farther than 10 Mpc 
at the 99.9% level. At this level of confidence, then, none of the bursts that 
have been observed by BATSE are in nearby galaxies. A signature from the 
Coma cluster or the "Great Wall" (~ 70 Mpc) is not likely to be seen in 
the data at present, since only a few bursts could have originated from these 
distances. 

The median distance to the nearest burst is 170 Mpc. Since the bright- 
est burst in the 3B catalog has a fluence of 7.8 x 10~^ erg cm~^ in 7-rays, 
this implies that the total energy released in 7-rays during a burst event is of 
order 3 x 10^^ h~'^ erg. 

As the number of observed 7-ray bursts keeps increasing, the distance limit 
will improve. In fact, with 3000 burst locations, the clustering of bursts 
might just be detectable and would provide compelling evidence for a 
cosmological origin. If it is not detected, the redshift to the "edge" of the 
burst distribution would be put at z ^ 1 or beyond. 
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